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Science, 
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NASA Marshall Space Flight Center 


Convection effect and source terms are two primary sources of difficulties 

in computing turbulent reacting flows typically encountered in propulsion 

devices. The present work intends to elucidate the individual as well as the 

collective roles of convection and source terms in the fluid flow equations 

and to devise appropriate treatments and implementations to improve our 

current capability of predicting such flows. A controlled variation scheme 

(CVS) has been under development in the context of a pressure-based 

algorithm, which has the characteristics of adaptively regulating the 

amount of numerical diffusivity, relative to central difference scheme, 

according to the variation in local flow field. Both the basic concepts and 

a pragmatic assessment will be presented to highlight the status of this 
work. 
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Overview of Algorithm 
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Treatment of Convection 

Convection terms: Strong Nonlinearity 
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Present Approach 


Q 

.ti H 

GO w 

O W) 

8 .5 S? 

£ -S > 

^ G O 

g s « 

•g Q I 

£ G -G 


5 £ 

S p2 

J2 W) 
uE G 


^ X5 
G U 


T3 X) 

a> *53 

IS oo 

“ <D 

5 g< 

co S 

C/3 r g 

<U U 

CD 

C3 

U A 


1893 


Harten’s Implicit TVD Scheme 
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Linear Steady Burgers’ Equation 
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Critical Value : I P I > 2 
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Linear Steady Burgers’ Equation With Source 
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Model Problem II: Different Schemes 
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Pcvs 
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Model Problem II: Dirichlet B.C., P=100 
Effective Cell Peclet Number 
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Treatment of Source Terms 
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TREATMENT OF SOURCE TERMS: A LONGITUDINAL 
COMBUSTION INSTABILITY PROBLEM 

Besides convection terms, source terms (if present) in the Navier-Stokes 
equations can be strong enough to cause numerical difficulties such as a loss 
of accuracy in the form of spurious oscillations in the solution profiles or 
numerical instability. This is so because strong source terms can be 
sufficiently stiff and the time and length scales imposed by them may not 
be commensurate with those imposed by convection, for example. Thus, 

due attention has to be paid to the source terms and not just to the convection 
terms. 

A one-dimensional longitudinal combustion instability problem is chosen 
which has a strong heat release source term. The high accuracy TVD type 
of convection treatment in a sequential solver (second fig. clockwise: top 
right) is seen to provide higher accuracy than the first-order upwind 
scheme (first figure: top left), as evident from the amplitudes of the ten 
pressure modeshapes shown in the viewgraph. However, the TVD type of 
convection treatment without any special source term treatment yields 
spurious oscillations in modehapes numbered 5, 6 and 7 (second figure 
clockwise). From the corresponding heat release modeshapes (third fig. 
clockwise), it is clear that modes 5, 6 and 7 are the modes of maximum heat 
release, thus demonstrating that when source terms become stiff enough 
they may lead to spurious oscillations. This can be resolved by increasing 
the amount of numerical dissipation in the scheme (by varying 6) but this 
is accompanied by an overall smearing of solution profiles. However, 
special source term treatment such as MacCormack’s predictor— corrector 
method or Strang’s time-splitting method (here, the latter) can resolve the 
problem by suppressing any spurious oscillations without the need of any 
extra numerical damping. This is clearly evident from the bottom left plot 
(fourth fig. clockwise). 


1901 


Special Source Term Treatment 


• Conservation law with a source term 

w t + f(w) x = ♦.("') 


• Treatment 

a) MacCormack’s Predictor-Corrector Method 

b) Operator Splitting (Strang’s Time-Splitting) 

W n+1 = Sy(At/2) S f (At) S f (At/2) W n 

where S f represents the numerical solution operator for 

W t * f( W )x = 0 

and Sj, is the numerical solution operator for the ODE 

w t = \J/(w) 
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1-D Combusting Flow in a Duct 

• To illustrate the effect of a strong source term (heat 
release) on numerical accurcay 

® To demonstrate the efficacy of special source term 
treatment for a strong source term 



Basic Equations 


3p + 3(p u) = Q 
dt dx 


dm d(mu) dp 1 ^ T ll 

dt dx dx Re^ dx 


3E d(Eu) _ d(pu) 1 [ 1 d dh, 9 , , 

dt ax ax at at + at Tn 


+ Q (x,t) 
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Combustion Instability Problem : Mode Shapes 



TVD Scheme; Source Treatment TVD Scheme; No Source Treatment 

= 0 
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(d) Controlled variation scheme; Xr/H-4.43 

Streamfunction: Backward-Step, Turbulent Flow 

Re = 132,000 
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TVD Schemes can be Effectively Extended to Sequential Solvers 


d) 

t3 

.a 

a 

o 

9 

g 

o 

*-G 

o 

d> 

> 

g 

0 

U 

V-i 

d> 

1 

£ 

33 

i— H 

O 

<D 

P-i 

13 

L) 

*G 

bi) 


S g 

O o 


LL, 
<D ^ 
> bD 

G G 
o \G 
G 


4-h 3 

w 8 

d) 

U Pt 


g 

£ 

Vx 

d> 

1 

z 

4-» 

d) 

i— ( 

o 

d) 

P< 

13 

U 

13 

o 

o 

P 

d> 

> 


G 

o 


O G 
d) Cd 

to ' 

w 

CO Q 
d) . 


GO 

C/D 


o 

G 

T3 

d> 

P< 

t 


G 

O 

• H 

<3 


G 

. 2 

rf ’’O 
.53 

§ o 

*G G 

s> G 

J£} o 

S o 


x) 

• rH 

Vx 

o 

jg 

4-> 

£ 

13 

o 

* 

G 

O 

+-> 

o 

G 

O 


CO 

d> 

CO 

G 

£ 

o 

d> 

Q 

G 

O 

• rH 
+-> 

G 

Px 

• rH 

C/D 

C/D 

• 

Q « 
13 § 

.a a 

U fly 

§ G 
S iG 

G <D 


G 

G 

£ 

£ 

r* 

bb 

d> 

r» 

CO 

’g 

Ph 

d> 

o 

G 

Vx «i-x 

d> '-zJ 
4-x ’G 
4-x fG 
•jr* G 

Q N 

-rH 

'S S 

G P6 

d) r-X 

u & 


o 

rH 

C/D 


G 

<D ^ 

> I 3 

^ P - ^ 

_j CO 
^3 d) 

erf O 

d> fa 
oo O 

I — > 4_( 

P G 

>> w 

r-H 

s, t 

“ c/d 
00 £ 

> ,2 

U Gh 


1908 



REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No. 0704-0188 

Public reporting burden for this collection of information n estimated to average t hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and competing and reviewing the collection of information Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for information Operations and Reports, 1215 Jefferson 
Davts Highway, Suite 1204, Arlington. VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 

1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

July 1993 Conference Publication 

4. TITLE AND SUBTITLE 

Eleventh Workshop for Computational Fluid Dynamic Applications in 
Rocket Propulsion — Part II 

5. FUNDING NUMBERS 

C. AUTHORS) 

R.W. Williams, Compiler 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 

B. PERFORMING ORGANIZATION 
REPORT NUMBER 

M-726 

9. SPONSORING /MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546 

10. SPONSORING /MONITORING 
AGENCY REPORT NUMBER 

NASA CP-3221 

fl. SUPPLEMENTARY NOTES 

Prepared by Structures and Dynamics Laboratory, Science and Engineering Directorate. 

12a. DISTRIBUTION /AVAILABILITY STATEMENT 

Subject Category: 34 
Unclassified — Unlimited 

12b. DISTRIBUTION CODE 

13. ABSTRACT (Maximum 200 words) 


Conference publication includes 79 abstracts and presentations and 3 invited presentations 
given at the Eleventh Workshop for Computational Fluid Dynamic Applications in Rocket Propulsion 
held at George C. Marshall Space Flight Center, April 20-22, 1993. The purpose of the workshop is to 
discuss experimental and computational fluid dynamic activities in rocket propulsion. The workshop is 
an open meeting for government, industry, and academia. A broad number of topics are discussed 
including computational fluid dynamic methodology, liquid and solid rocket propulsion, 
turbomachinery, combustion, heat transfer, and grid generation. 


14. subject terms S p ra y^ injector, computational fluid dynamics, rocket propul- 
sion, liquid rocket, solid rocket, turbopump, turbomachinery, combustion, 
methodology, impeller, inducer, heat transfer, grid generation, nozzle, plume 

15. NUMBER OF PAGES 

955 

16. PRICE COOE 

A99 

17. SECURITY CLASSIFICATION 
OF REPORT 

18. SECURITY CLASSIFICATION 
OF THIS PAGE 

19. SECURITY CLASSIFICATION 
OF ABSTRACT 

20. LIMITATION OF ABSTRACT 

Unclassified 

Unclassified 

Unclassified 

Unlimited 


NSN 7540-01-280-5500 


Standard Form 298 {Rev. 2-89) 

Prescribed by ANSI Std Z39-18 
298-102 



















